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Abstract 

We report the proof that the extension of Gibrat's law in the middle scale region is unique 
and the probability distribution function (pdf) is also uniquely derived from the extended 
Gibrat's law and the law of detailed balance. In the proof, two approximations are employed. 
The pdf of growth rate is described as tent-shaped exponential functions and the value of 
the origin of the growth rate distribution is constant. These approximations are confirmed 
in profits data of Japanese companies 2003 and 2004. The resultant profits pdf fits with the 
empirical data with high accuracy. This guarantees the validity of the approximations. 
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In the large scale region of wealth, assets, sales, profits, income, the number of employees 
and etc (x), the cumulative probability distribution function (pdf) P(> x) obeys a power-law 

P(> x) oc for x > xo , (1) 

where xq is a certain threshold. This power-law is called Pareto's law [1], which has been well 
investigated by using various models. 

A few years ago, Fujiwara et al. [2] show that Pareto's law is derived from Gibrat's law [3] 
and the law of detailed balance by using no model. Gibrat's law is observed in the large scale 
region. It argues that the conditional pdf of growth rate Q(R\x±) is independent of the initial 
value x\\ 

Q{R\ Xl ) = Q{R) . (2) 

Here growth rate R is defined as R = xijx\ and Q{R\x\) is defined by using the pdf P{x\) and 
the joint pdf Pir{x\,R) as Q{R\x\) = P\r{x\, R)/P(xi). The detailed balance observed in a 
stable economy is time-reversal symmetry (x± <-> X2): 

Pi 2 (xi,x 2 ) = Pi2(x 2 ,xi) . (3) 

Here x\ and x 2 are two successive wealth, assets, sales, profits, income, the number of employees, 
etc. and Pv2.(x\,x<2) is the joint pdf. 
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This derivation is significant for the understanding of the mechanism of Pareto's law, how- 
ever, it is valid only in the large scale region where Gibrat's law (2) holds. It is well known 
that Pareto's law is not observed below the threshold xq [3, 4]. The reason is thought to be the 
breakdown of Gibrat's law [2], [5]-[7]. 

In Ref. [8], Gibrat's law is extended in the middle scale region by employing profits data of 
Japanese companies in 2002 and 2003. We approximate the conditional pdf of profits growth 
rate as so-called tent-shaped exponential functions 

Q(R\ Xl ) = d( Xl ) R-t+i^)- 1 for R > 1 , (4) 
Q(R\ Xl ) = d( Xl ) fl+M*i)-i for R < 1 . (5) 

By measuring t± we have assumed the x\ dependence to be 

t±(x\) = t±(x ) ±a± In — , (6) 

and have estimated the parameters as [9] 

a + ~ a_ ~ for x\ > xq , (7) 

a + ~ a_ 7^ for x m j n < x\ < xq , (8) 

t+(x ) -t-(x ) ~ ii . (9) 

From extended Gibrat's law (6) - (9) and the detailed balance (3), we have derived the pdf in 
the large and middle scale region uniformly as follows 

P{x) = Cx~^ +1) e' aln2 ^ for x > x min , (10) 

where a = (a + + a_) /2. This is confirmed in the empirical data. 

In this study, we prove that the x\ dependence of t± (6) with a + = a_ is unique if the 
pdf of growth rate is approximated by tent-shaped exponential functions (4), (5). This means, 
consequently, that the pdf in the large and middle scale region (10) is also unique if the x\ 
dependence of d{x\) is negligible. We confirm these approximations in profits data of Japanese 
companies 2003 (xi) and 2004 (X2) [10] and show that the pdf (10) fits with empirical data 
nicely. 

In the database, Pareto's law (1) is observed in the large scale region whereas it fails in 
the middle one (Fig. 1). At the same time, it is confirmed that the detailed balance (3) holds 
(Fig. 2). The scatter plot in Fig. 2 is different from one in Ref [8]. The reason is that the 
identification of profits in 2002 and 2003 in Ref. [8] was partly failed. As a result, the pdfs of 
profits growth rate are slightly different from those in this paper. The conclusion in Ref. [8] is, 
however, not changed. 

The breakdown of Pareto's law is thought to be caused by the breakdown of Gibrat's 
law in the middle scale region. We examine, therefore, the pdf of profits growth rate in the 
database. In the analysis, we divide the range of x\ into logarithmically equal bins as xi £ 4 x 
j 10 i+o.2(n-i) ^ 1Q i+o.2n] thousand yen with n = 1, 2, • ■ ■ , 20. In Fig. 3, the probability densities for 
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r are expressed in the case of n = 1, • • • , 5, n = 6, • ■ ■ , 10, n = 11, • ■ ■ , 15 and n = 16, • ■ ■ , 20, re- 
spectively. The number of the companies in Fig. 3 is "22,005", "89,507", "85,020" and "24,203", 
respectively. Here we use the log profits growth rate r = log 10 R- The probability density for r 
defined by q(r\x\) is related to that for R by log 10 Q{R\x\) + r + log 10 (lnl0) = log 10 g(r|xi) . 
From Fig. 3, log 10 q(r\x\) is approximated by linear functions of r as follows 



logio^^i) 
\og w q{r\xi) 



c(xi) - t + (xi) r 
c(xi) + t-(xi) r 



for r > , 
for r < . 



(11) 
(12) 



These are expressed as tent-shaped exponential functions (4), (5) by d(x\) = lO^^/ln 10 • In 
addition, the x\ dependence of c{x\) is negligible for n = 9, • • • , 20 (Fig. 3). The validity of these 
approximations should be checked against the results. 

We show that the x\ dependence of t± (6) is unique under approximations (4), (5). By the 
use of the relation of Pi 2 (xi, x 2 )dx\dx 2 = P\r{x\, R)dx\dR, the detailed balance (3) is rewritten 
as P\r{x\,R) = R^ 1 Pir(x2, R^ 1 ) ■ Substituting the joint pdf P\r(x\,R) for the conditional 
probability Q(R\x\), the detailed balance is expressed as 

P( Xl ) 1 Q(i?- 1 |x 2 ) 



P(x 2 ) R Q{R\xi) ' 
Under approximations (4) and (5), the detailed balance is reduced to 



P(xi) _ d(x 2 ) 



U+t+(xi)-t_(x 2 )+l 



(13) 



(14) 



P(x 2 ) d(xi) 

for R > 1. By expanding Eq. (14) around R = l, the following differential equation is obtained 



l + t+(x) 



P(x) + x P (x) = 0, 



(15) 



where x denotes X\ and P(x) = P{x)d(x). The same differential equation is obtained for 

where t + (x) — t-{x) = g{x) and 



Cx~ 



-G(x) 



R < 1. The solution is given by P(x) 
J g(x)/x dx = G(x). 

In order to make the solution around R = 1 satisfies Eq. (14), the following equation must 
be valid for all R: 



-G(x) + G(Rx)= t+(x)-t-(Rx) 



InR . 



(16) 



By expanding the derivative of Eq. (16) with respect to x around R = 1, following differential 
equations are obtained 



t+ (x) + t- (x) + t+ (x) + 1- (x) = , 



2 t+ (x) + t. 
The solutions are given by 

t+(x) = - 



3x t. 



X 



t + W(x) + 2t. 



(3), 



. 



C_2 , _> 



In" x + (C+i - C_i) lnx + (C +0 - C_ ) , 



C 

t_( x ) = — ^ln 2 x + C_ilna; + C_o • 



(17) 
(18) 

(19) 
(20) 
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To make these solutions satisfy Eq. (16), the coefficients must be C_2 = and C+i = 0. Finally 
we conclude that t±(x) is uniquely expressed as Eq. (6) with a + = a_. 
Under approximations (4) and (5), we obtain the profits pdf 

P{x) = P(x)d{x) = Cx~ { ^ e~ ahl2 ~h , (21) 

where we use the relation (9) confirmed in Ref. [9]. In Fig. 4, t± hardly responds to x\ for 
n = 17, • • • , 20. This means that Gibrat's law holds in the large profits region. On the other hand, 
t+ linearly increases and t_ linearly decreases symmetrically with log 10 x\ for n = 9, 10, • ■ ■ , 13. 
The parameters are estimated as Eq. (7) and (8) with a (= q + = a_) ~ 0.14 and xq = 
4 x 10 1+0 - 2 ( 17_1 ) ~ 63,000 thousand yen. Because the x\ dependence of c{x\) is negligible in 
this region, the profits pdf is reduced to Eq. (10). We observe that this pdf fits with the empirical 
data nicely in Fig. 5. Notice that the estimation of a in Fig. 4 is significant. If we take a slightly 
different a, the pdf (10) cannot fit with the empirical data (a = 0.10 or a = 0.20 in Fig. 5 for 
instance). 

In this paper, we have shown the proof that the expression of extended Gibrat's law is 
unique and the pdf in the large and middle scale region is also uniquely derived from the 
extended Gibrat's law and the law of detailed balance. In the proof, we have employed two 
approximations that the pdf of growth rate is described as tent-shaped exponential functions 
and that the value of the origin of the growth rate distribution is constant. These approximations 
have been confirmed in profits data of Japanese companies 2003 and 2004. The resultant pdf 
of profits has fitted with the empirical data with high accuracy. This guarantees the validity of 
the approximations. 

For profits data we have used, the distribution is power in the large scale region and log- 
normal type in the middle one. This does not claim that all the distributions in the middle 
scale region are log-normal types. For instance, the pdf of personal income growth rate or sales 
of company is different from tent-shaped exponential functions [2]. In this case, the extended 
Gibrat's law takes a different form. In addition, we describe no pdf in the small scale region [11]. 
Because the X\ dependence of d{x\) in the small scale region is not negligible (Fig. 3). Against 
these restrictions, the proof and the method in this paper is significant for the investigation of 
distributions in the middle and small scale region. We will report the study about these issues 
in the near future. 
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Figure 1: Cumulative probability distributions Figure 2: The scatter plot of all companies in 

P(> x\) and P(> X2) for companies, the profits the database, the profits of which in 2003 (xi) 

of which in 2003 (x±) and 2004 (#2) exceeded and 2004 (#2) exceeded 0, x\ > and x^ > 0. 

0, x\ > and X2 > 0. The number of the companies is "227,132". 
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TheLogarithmofProfitsGrowthRater=log 10 x ,/x , 



-2-1012 
TheLogarithmofProfitsGrowthRater=log 10 x ,/x , 




-2-1012 -2-1012 

TheLogarithmofProfitsGrowthRater=log 10 x ,/x , TheLogarithmofProfitsGrowthRater=log 10 x ,/x , 



Figure 3: The probability density distribution q(r\x\) of the log profits growth rate r = 
log 10 X2/X1 from 2003 to 2004. The data points are classified into twenty bins of the initial prof- 
its with equal magnitude in logarithmic scale, x\ £ 4 x [I0 1+a2<n-1 ), io 1+0 - 2ri ] (n = 1, 2, ■ ■ ■ , 20) 
thousand yen. 
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Figure 4: The relation between the lower bound 
of each bin x\ and t±(x±). From the left, each 
data point represents n = 1, 2, • • • , 20. 



Figure 5: The probability distribution function 
(pdf) P{x2) for companies, the profits of which 
in 2003 (ari) and 2004 (x 2 ) exceeded 0. 



